Adiabatic noise-induced escape rate for nonequilibrium open systems 
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We consider the motion of an overdamped particle in a force field in presence of an external, 
adiabatic noise, without the restriction that the noise process is Gaussian or the stochastic process 
is Markovian. We examine the condition for attainment of steady state for this nonequilibrium open 
system and calculate the adiabatic noise-induced rate of escape of the particle over a barrier. 
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I. INTRODUCTION 

The motion of a Brownian particle in a fluid was first 
explained by Einstein in terms of fast thermal motion of 
fluid molecules striking the Brownian particle and caus- 
ing it to undergo a random walk. One essential require- 
ment of the theory is that the noise is of internal ori- 
gin. This implies that the dissipative force which the 
Brownian particle experiences in course of its motion in 
the fluid and the stochastic force acting on the particle 
as a result of random impact of molecules arise from a 
common mechanism. From a microscopic point of view 
the system-reservoir Hamiltonian description developed 
over the last few decades ||,^ suggests that the cou- 
pling of the system and the reservoir co-ordinates de- 
termines both the noise and the dissipative terms of the 
Langevin equation describing the motion of the particle. 
It is therefore not difficult to anticipate that these two 
entities get related through the celebrated fluctuation- 
dissipation theorem (Einstein's relation for diffusion and 
mobility is the first of its kind). The spiritual root of 
fluctuation-dissipation relation lies at the dynamical bal- 
ance between inward flow of energy due to fluctuation 
of the reservoir into the system and the outward flow of 
energy from the system to the reservoir due to dissipa- 
tion of the system. This dynamical balance is essential 
for attainment of a steady state of the nonequilibrium 
system. In their treatise on nonequilibrium statistical 
mechanics Lindenberg and West |^ have classified these 
systems as thermodynamically closed. However, there 
are quite a large number of physical situations, a com- 
prehensive account of which has been given in , where 
the system is thermodynamically open, i.e., when the 
system is driven by an external noise which is indepen- 
dent of system's characteristic damping. These are im- 
portant for describing the effects of pump fluctuations on 
the emission of a dye laser , effects of ffuctuating rate 
constants on a chemical reaction , effects of noise on 
electronic parametric oscillators etc. The dynamics is 
still governed by a Langevin equation. Although there 
exists no relationship between the ffuctuation and the 
dissipation in such a situation it is interesting to search 
for the condition under which the steady state of the 



thermodynamically open system is attained. Physically 
the sytem is described by three timescales ; the timescale 
of dissipation or relaxation /3~^ unlike the closed system 
is independent of correlation time Tc of the system. In 
the present communication we consider the system to 
be driven by fluctuation which is adiabatically slow such 
that it is characterized by a very long correlation time Tc 
where 



/3-^ < At < Tc 



(1) 



Here At refers to the timescale over which we look for 
the average motion of the system. The latter inequality 
implies that l3~^, i.e., the inverse of damping constant 
defines the shortest timescale in the dynamics in con- 
trast to the case of standard Brownian dynamics which 
obeys At Tc. Secondly, we do not put the 

restriction that the stochastic processes is Markovian or 
the noise is Gaussian. The two assumptions have been 
discussed so much in the recent literaure that it is neces- 
sary to emphasize that the present consideration is free 
from these assumptions. We mention, in passing, that 
some interesting limitting situations have been examined 
for linearized systems driven by Gaussian noise processes 



Our aim in the present article is two-fold. First, we ex- 
amine the condition for a steady state distribution of the 
system driven by an external noise obeying the inequality 
(0) . Second, Having obtained this condition we calculate 
the external, adiabatic noise induced steady state escape 
rate over a barrier in the spirit of Kramers'-Smoluchowski 
theory. We show how the third-order noise plays an im- 
portant role in both of these issues of stochastic dynam- 
ics. 



II. THIRD ORDER NOISE AND STEADY STATE 
PROBABILITY DENSITY FOR OPEN SYSTEMS 

To start with we consider the equation of motion of a 
particle of unit mass in one dimension when it is acted 
upon by an external field of force corresponding to a po- 
tential V{x) and an external, adiabatic stochastic force 
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where p and the correlation time Tc of ^(t) satisfy the 
inequality ([|) . Also note that by virtue of this we have 
considered the ovcrdamped limit. In a preceding paper 
1^ the equation of motion for probability density distri- 
bution function P{x, t) in phase space corresponding to 
the Langevin description (|) was derived. It has been 
shown that P{x,t) obeys the differential equation of mo- 
tion which contains third order terms (beyond the usual 
Fokker-Planck terms) giving rise to third order noise. 
The appearance of these terms is generic for the stochas- 
tic process pertaining to the separation of timescales (Q) 
we consider here. The general expression for time evolu- 
tion of probability density function is given by B 



dt 



dx^ 



Co, Ci and C2 in Eq.(^ are as follows : 
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where we have also assumed, for convenience {^{t)) = 0. 

The equation describes the time evolution of an 
overdamped particle in a force field (derivable from a 
potential V{x)) simultaneously subjected to an external 
adiabatic stochastic force, co, ci and C2 measure the 
strength of the noise term. While the first term in Eq.(|^) 
can be identified as the usual deterministic dynamical 
term, the second and the third terms involve second 
and third order diffusion coefficients due to stochastic- 
ity of ^(t). The expansion scheme associated with Eq.(||) 
should not be confused with Kramers-Moyal [KM] ex- 
pansion (consideration of third order diffusion term of a 
KM expansion may lead to serious interpretive difhcul- 
ties Q because the probability distribution function often 
turns out to be negative) which serves as the standard 
starting point for analysis of stochastic processes with 
fast noise corresponding to the separation of timescale 
/3~^ ^ Ai ^ Tc. What is implicit in a KM expansion 
is that the moments are assumed to be linear in r and 
thus the validity of the coefficients of a KM expansion 
rests on the smallness of r. The cummulant expansion 
similarly relies on the smallness of correlation time 
Tc- On the other hand Eq.^ is based on "adibatic fol- 
lowing approximation" which involves an expan- 
sion in 1//3 which is evident from the nature of 
the coefficients of the terms in the right hand side. The 
remarkable departure from the standard form of Fokker- 
Planck equation (Smoluchowski equation) is due to the 



presence of the third order noise. For other details we 
refer to 

In the next step we recast the third order equation 
(^) in the form of the familiar continuity equation and 
identify a steady state current J in the following equation 
for a steady state probability distribution function P{x) 
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where a and b are given by 
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We consider a Kramers type potential V{x) shown in 
Fig.(l) and look for the steady state current at the bar- 
rier top by linearizing the potential at a; = 0. If ujq refers 
to frequency of the inverted well and -Bo defines the po- 
tential at the barrier top, one obtains the following inho- 
mogenous modified Bessel equation of order v 
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In the above derivation following transformations and ab- 
breviations have been used 



P{x) = x^^-'''>/^ W{x) , C = V6a 
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The homogenous counterpart corresponding to the above 
Eq. (^ is the standard modified Bessel equation of order 
v. The general solution of Eq.(|^) can be written as 



WiC) = AUO+BK,iC) 
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where and i^i/(C) are modified Bessel functions of 

order v; A and B are the two arbitrary constants of inte- 
gration corresponding to the homogenous part of Eq.(0). 
The D containing term results from the particular inte- 
gral corresponding to the inhomogenous contribution of 
Eq.(^ obtained by the method of variation of parame- 
ters. Making use of the relations (|8^) we revert back to 
the original variables x and P{x) to obtain the general 
solution of Eq. (^ as 

P{x) = A x^"" luiVbx) + B x-" K^Vbx) 
h{Vbx) / x^K 
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y/bx 



x'''U{Vbx')dx' 



(10) 
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We now impose two natural boundary conditions on 
the solution for the probability density function ; (i) P{x) 
vanishes for |x| — > oo and (ii) P{x) remains finite for all 
X. The first condition ascribes the essential notion of a 
probability function which demands A = 0. The second 
one is a necessary requirement since the relevant modified 
Bessel function is singular at a; = 0. This singularity at 
a; = can be removed by constructing a relation between 
the constants B and D of the form 



B= (-l)"i?-^2(2"+i)/2 
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such that all powers of j vanishes identically from the 
exact solution of the modified Bessel equation (|^) 
which satisfy both the boundary conditions simultane- 
ously provided h' is an odd half integer, i.e., 

i^ = n+ (1/2) , 71= 1,2,... 



so that from (g) and the definition of in (8b) we have 
(note that n = is not physically allowed) 
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The above Eq.(|l2|) relates the dissipation constant /3 to 
the correlation of fluctuations ( coi , C2 ) of the external 
adiabatic noise as defined in Eq. (4) . This is a fluctuation- 
dissipation like relation for the open systems. The in- 
tegers n characterize the distinct stable steady states 
in terms of the positive definite probability distribution 
function as given by (for a; > 0) 



P„+i(x) = (-l)"D2" n\ 
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The normalization constant D is related to the steady 
state current J (positive, implying a flow from left to 
right) through the first of the relations (|^) and is given 
by 



D = b 
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where A is large but finite and the function Ein(x) is 
defined llSl as 



Ein(x)=^(-1)'= 
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For X < 0, the corresponding probability distribution 
function can be calculated from the symmetry of the dif- 
ferential equation (|^). 

Expression ( p^ ) reduces to the following simple form 
for, say, n = 1, which is depicted in Fig. (2) for several 
values of third order noise strength C2 
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Thus the steady states are physically allowed only when 
one can relate the dissipation (3 to the strength of sec- 
ond and third order external noise for these integers n. 
Our analysis suggests that such a steady state condition 
( p^ ) for the open systems can be realized at least in the 
specific issues as in the present instance. The third order 
noise is an essential ingredient for the dynamic balance 
we have referred to. It is thus also apparent why the ther- 
modynamic open systems described by a Fokker-Planck 
equation for fast fluctuations which include only second 
order diffusion coefficient reaches a steady state with no 
fluctuation-dissipation relation. 



III. DYNAMICS OF BARRIER CROSSING 
INDUCED BY ADIABATIC NOISE 

Once the condition for attainment of the steady state is 
realized it becomes possible to consider the situation such 
that a particle in a force field, e.g., originally confined in 
a potential well may escape under the influence of exter- 
nal adiabatic noise by maintaining a steady state current 
over the barrier. It is therefore pertinent to calculate the 
noise-induced escape rate in the spirit of Kramers and 
Smoluchowski and to elucidate the aspect of dependence 
of escape rate on dissipation. The counterpart of the lat- 
ter issue in the theory of fast fluctuation is the wellknown 
turn-over problem in Kramers' theory |l^Jl7|] . 

To proceed further we again make use of a Kramers' 
type potential V{x) ( Fig.(l) ) under which the particle 
moves, in the third order equation (0) for the probability 
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distribution function. The popular flux-over-population 
method [ p^JT?! is then employed. The calculation rests 
on the evaluation of two quantities; (i) the steady state 
current J over the barrier top, (located at x = 0) that 
results if the particles are continuously fed into the do- 
main of attraction (say, in the region of left well) and are 
subsequently and continuously removed in the neighbor- 
ing domain of attraction, (ii) steady state population Ua 
in the initial domain of attraction, i.e., the left well. The 
rate is defined by 



K. = J/ua ■ 



(18) 



For linearization the potential V{x) at the bottom of 
the left well at a: = — A we approximate V{x) ~ ^ oj'^ {x+ 
A)^, where ojb refers to the frequency at the bottom of the 
left well. The above approximation to left well and J = 
condition reduce the third order equation of motion 
(Eq.(||)) for probability density Pb{x) near the bottom of 
the well to the following form. 

The above equation is valid near the bottom of the left 
well {x ~ — A). Here 7' is defined as 
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and a by F.q.@. Eq.(|T|) can be written as 



making use of the following sets of of transformations. 



y{z) = W{z) , ly' = 
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The solutions of Eq.(|T]) are again the modified Bessel 
functions. Reverting back to original variables, the gen- 
eral solution for the steady state probability distribution 
near the bottom of the left well is given by 

Pb[x) ^A'{x + ^Y' I,, [Vb (x + A)] 

+B' {x + Ay' K,, [Vh {x + A)] , (23) 

A' and B' are the two arbitrary constants of integration. 
The solution Pb{x) must satisfy the following boundary 
conditions ; (i) Pb{x) must vanish at x — > cxd and (ii) 

n+- 

Pb{—A) = ^ (^^)i where the stationary probability 

n+- 

Pf. ^ (—A) corresponds to the vanishing current J ~ 
along X pertaining to the homogenous version of Eq. (^ . 

As usual, Pf. ^ {x) must also satisfy the boundary con- 

n+- 

dition that for |a;| ^ 00, ^ {x) vanishes. The first 



condition leads to A' = and the second one gives the 
value of B' in terms of B of Eq. ([lO|) and therefore 



Pb{x) ^B' {x + ay' K,, [V6 {x + A)] 



(24) 



where 
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Therefore Pb{x) in Eq.(p4]) may be expressed as 
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The above distribution which is valid near the bottom 
of the left well may be used to calculate the population 
inside the left well as. 



.0 

la = 2 Pb{x) dx . 



(27) 



Due to the presence of Ky{x) the probability Pb{x) is a 
rapidly decreasing function. We may extend the above 
integration limit to infinity. This yields ( using Eq.(24) ) 
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Using the relations (^5|) and ( 



we finally have 
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Having determined the population of the left well, Ua 
and the steady state current, J ( from first of the relations 
m (1^) and (|l|) ) over the barrier we are now in a position 
to calculate the escape rate /C„+i (= J I no). We quote 
the final result : 
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A-s refer to the zero's of the potential V{x) as shown in 
Fig.(l) and by virtue of hnearization of V{x) at x = ; 
A is approximately given by A = (2£'o/wo)^^^ ■ 

The above expression can be made more transparent 
by demonstrating a representative transition rate, say, 
for n = 1 as follows : 
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The above expression is analogous to Kramers' formula 
for the rate of escape from a potential well over a finite 
barrier of height Eq under the influence of an external 
nonthermal adiabatic noise pertaining to the timescale 

The escape rate expressions derived above suggest 
that the rate approaches zero both for /3 large and 
(3 small. This behavior is somewhat reminiscent of 
Kramers' theory ||l^, where it was noted earlier that 
this two limiting behavior implies a maximal rate at 
some damping value (3. The rate therefore undergoes 
a turnover in a form of a bell-shaped curve. In Fig. (3) 
we plot a representative variation of the escape rate ver- 
sus dissipation (3 for different third order noise strength. 
With increasing friction, the rate undergoes a turnover 
from an increasing behavior at low friction to an inverse 
behavior in the high friction limit. 



IV. CONCLUSION 

In conclusion, we consider the motion of particle in 
a force field, simultaneously subjected to adiabatic fluc- 
tuations of external origin. The equation of motion for 
probability distribution function includes a third order 
noise term. We show that although the system is ther- 
modynamically open, the specific interplay of the char- 
acteristic dissipation of the system and the correlation of 
fluctuations due to external non-Gaussian noise leads to 
distinct steady states for the open system. We calculate 
the external adiabatic noise-induced rate of escape of a 
particle confined in a well. A typical variation of the es- 
cape rate as a function of dissipation which is reminiscent 
of Kramers' turn-over problem, has been demonstrated. 
In view of several experimental investigations on exter- 
nal noise induced transitions in the recent past [||-^ , the 
study of thermodynamically open systems has been spe- 
cially relevant in both physical and chemical sciences. Al- 
though the driving noise processes in these cases are fast, 
suitable extension to adiabatic noise limit may throw new 
light on the present issue. The population inversion in 
a two-level system by an adiabatically varying stochastic 
electromagnetic field, as considered earlier in Refs. [ p^Jl9| 
may also be worthwhile candidate for further studies in 
this context. 
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Figure Captions 

1. Fig.(l) : A schematic plot of the Kramers' type 
potential V{x) used for the calculation of barrier 
crossing dynamics. 

2. Fig. (2) : The normalized probability distribution 
function ^3/2(2;) is plotted as a function of a; for 
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various values of the third order noise strength C2 
( /3= 1.0 and A = 2.5 ). 

3. Fig. (3) : Escape rate /C3/2 is plotted as a function 
of the characteristic dissipation (3 of the system for 
various values of C2 ( cqi = 7.0, ujb = 0.80 and 
A = 2.5 ). 



6 



V(x) 




